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! Abstract 

• We calculate the cross sections for the neutral Higgs-boson production at 

D , e^e^ colliders in the Minimal Supersymmetric Standard Model (MSSM) using the 

Feynman-diagrammatic approach and the on-shell renormalization scheme. We in- 
corporate the Higgs-boson propagator corrections, evaluated up to two-loop order, 
^ . into the prediction of the cross sections for the Higgs-boson production mechanism 

^ I e~^e~ — > hZ, hA. The propagator corrections consist of the full one-loop contribu- 

tion, including the effects of non-vanishing external momentum, and at the two-loop 
level of the dominant corrections of 0{aas) and further sub-dominant contributions. 
The results are supplemented with the complete set of one-loop vertex and box cor- 
rections. The effects of the two- loop propagator corrections are investigated in 
detail. We briefly discuss also the effect of the box contributions for high ^/s. We 
compare our results with the case where only the corrections to the effective mix- 
ing angle, evaluated within the renormalization-group-improved one-loop Effective 
Potential approach, are taken into account. We find agreement better than 10% for 
LEP2 energies and deviations larger than 20% for = 500 GeV. 
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1 Introduction 



The search for the hght neutral Higgs boson is a crucial test of Supersymmetry that can 
be performed with the present and the next generation of accelerators. The prediction of 
a relatively light Higgs boson is common to all supersymmetric models whose couplings 
remain in the perturbative regime up to a very high energy scale Finding the Higgs 
boson is thus one of the main goals of high-energy physics. Possible indications for a Higgs 
boson with a mass around 115 GeV have recently been observed in the LEP analysis 0. 
As a second step, once a scalar particle has been found, it is necessary to measure the 
dominant production cross sections as well as the decay widths and branching ratios of 
the main decay channels to a high accuracy. Such precise measurements of Higgs-gauge- 
boson and Higgs- Yukawa couplings allow to experimentally investigate the details of the 
Higgs mechanism. Finally, also the Higgs self-couplings will have to be measured in order 
to reconstruct the Higgs potential. A future linear e^e~ collider with high luminosity 
can provide a sufficiently clean environment to measure both the Higgs couplings to 
other particles as well as its self couplings [§] with high precision, allowing to distinguish 
between a standard and a non-standard Higgs boson. In this paper we concentrate on the 
production of the neutral Higgs bosons of the Minimal Supersymmetric Standard Model 
(MSSM) in e+e~ annihilations and provide results for the production cross sections for 
general model parameters, together with a detailed numerical study. 

In the MSSM, the mass of the lightest Higgs boson, M^, is bounded from above by 
Mh ^ 135 GeV, including radiative corrections up to two-loop order p|-p^. The most 
promising channels for the production of the supersymmetric neutral Higgs particles at 
LEP2 energies and in the first phase of future e~^e~ colliders are the Higgs-strahlung 
process [|14 

e+e-^Zh{H), (1) 
and the associated production of a scalar and a pseudoscalar Higgs boson 

e+e-^Ah{H). (2) 

We compute the MSSM predictions for the cross sections of both channels in the 
Feynman-diagrammatic (FD) approach using the on-shell renormalization scheme. We 
take into account the complete set of one-loop contributions, thereby keeping the full 
dependence on all kinematical variables. The one-loop contributions consist of the correc- 
tions to the Higgs- and gauge-boson propagators, where the former contain the dominant 
electroweak one-loop corrections of OlGprnj), and of the contributions to the 3-point 



and 4-point vertex functions [|I^-|T9|. We combine the complete one-loop result with the 
dominant two-loop QCD corrections of O^GpcHsm^) |]12|,|13| and further sub-dominant 
corrections. In this way the currently most accurate results for the cross sections are 
obtained. 

Furthermore we show analytically that the Higgs-boson propagator corrections with 
neglected momentum dependence can be absorbed into the tree-level coupling using the 
effective mixing angle from the neutral CP-even Higgs-boson sector. We compare our 
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results for the cross sections with the approximation in which only the corrections to 
the effective mixing angle, evaluated within the renormalization-group-improved one-loop 
Effective Potential approach, are taken into account. For most parts of the MSSM pa- 
rameter space we find agreement better than 10% for the highest LEP energies, while for 
= 500 GeV the difference can reach 25%. 

The paper is organized as follows: in section ^ the basic formulae are presented and it 
is shown analytically how the Higgs-boson-propagator corrections with neglected external 
momenta are related to the effective mixing angle in the Effective Potential approach 
(EPA). Section ^ contains the numerical analysis for the production cross sections and 
the comparison of our full result with the EPA. The conclusions can be found in section ^. 



2 Cross sections for Higgs-particle production in e^e 
collisions 

2.1 Classification of radiative corrections 

The two Higgs-field doublets giving rise to electroweak symmetry breaking within the 



MSSM accommodate five physical Higgs bosons [0. At the tree-level, two input pa- 
rameters (besides the parameters of the Standard Model (SM) gauge-sector) are needed 
to describe the Higgs sector. We choose them to be tan/3, the ratio of the two vacuum 
expectation values, and Ma, the mass of the CP-odd Higgs boson. The CP-even neutral 
mass eigenstates are obtained from the interaction eigenstates by the rotation 





cos a sm a 
— sin a cos a 




(3) 



with the tree-level mixing angle a related to tan (3, Ma and Mz by 



tan2a = tan2/3 ^^_^f , < « < 0. (4) 



Two main sources for the production of supersymmetric Higgs particles in e+e colli- 
sions are the Higgs-strahlung process 



e+e" ^ ZHi, i = l,2 (5) 

(using the compact notation of eq. (|^)) and the associated production of scalar and pseu- 
doscalar Higgs bosons, 

e+e~ ^ AHi, z = 1,2. (6) 

We do not discuss here the possibility of Higgs-particle production by bremsstrahlung off 
heavy quarks (e.g. e+e" — bbHi, which can be significant for large tan/? [^), or by 
W^W~ fusion, which becomes important for center of mass system (CMS) energies of 
C»(500 GeV) 
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The set of diagrams taken into account for Higgs-strahlung e"'"e~ hZ is schematically 
shown in Fig. |l|, where a) is the tree-level diagram. The shaded blobs summarize the loops 
with all possible virtual particles, except photons in the Ze~^e~ vertex corrections^. More 
details can be found in Refs. |]16|,0,|19|. An analogous set has been evaluated for the 



second process, e~^e~ — > hA. 

For completeness, in Fig. ^ also contributions are shown (e.g. the A, G-Z mixing con- 
tributions and the longitudinal parts of the Z and 7-Z self-energies) that are proportional 
to the electron mass or vanish completely after contraction with the polarization vector 
of the Z boson. The different types of corrections can be summarized as follows: 

(i) Corrections to the e, Z, 7 and 7-Z self-energies on the internal and external lines 
and to the (initial state) Ze~^e~ and ■ye~^e~ vertices, b) - g). 

(ii) Corrections to the scalar and pseudoscalar propagators, h). 

(iii) Corrections to the ZZHi [ZAHi) vertex, i). 

(iv) Box-diagram contributions and t-channel-exchange diagrams, j) - 1). 
The corrections (i)-(iv) have a different relative impact: 

- Electroweak corrections of type (i) are typically of the order of a few percent (like 
in the Standard Model) and do not exhibit a strong dependence on any SUSY 
parameters. 

- The main source of differences between the tree-level and higher-order results are the 
corrections to the Higgs-boson self-energies (ii). They are responsible for changes 
in the physical masses and Mh and the effective mixing angle (via contri- 
butions to the renormalization constants, Z'^^^, for the external Higgs particles in 
the S'-matrix elements, see Sect. |2.2| ) predicted for given values of tan/5 and Ma- 
At the one-loop level these propagator corrections constitute the only source for the 
large correction of OlGpfn'f). At the two-loop level they exclusively give rise to 
contributions of OlGpasmf) and of (9(G|'mf). In this sense the propagator correc- 
tions define a closed subset of diagrams, being responsible for a numerically large 
contribution. 

- Corrections to the final-state vertices (iii) are typically larger than those of type (i), 
but smaller than the Higgs-boson propagator corrections. At LEP2 energies they 



can reach at most 7-10% [jT6[ for very low or very large values of tan/3, when the 



Yukawa couplings of the top or bottom quarks become strong. 

Finally, the box-diagram contributions (iv) depend strongly on the center-of-mass 
energy. They are of the order of 2-3% at LEP2 energies and may reach 20% for 
v/i = 500 GeV pf. 



^ These virtual IR-divergent photonic corrections constitute, together with real-photon 
bremsstrahlung, the initial-state QED corrections, which are conventionally treated separately and are 
the same as for the SM Higgs-boson production. 
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It should be noted that initial-state QED corrections as well as finite-width effects (allow- 
ing for off-shell decays of the Higgs and the Z boson) are not included in our calculation. 
However, by incorporating our result into existing codes, e.g. HZHA [Q, QED corrections 
and finite-width effects can automatically be taken into account]^. 



2.2 Higgs-boson masses and wave function renormalizations 

Radiative corrections induce mixing between the CP-even Higgs bosons, even if their 
mass matrix has been diagonalized at the tree levelQ. In the FD approach the higher- 
order corrected Higgs-boson masses, denoted by M/j, Af^, are derived by finding the poles 
of the /i, if-propagator matrix, which is given by 




(7) 



Here the Tj{q^) denote the renormalized Higgs-boson self-energies (throughout the pa- 
per we denote tree-level mass parameters by small letters and physical masses by capital 
letters). For these self-energies we take the result of the complete one-loop on-shell calcu- 
lation of Refs. |1^,|T6| together with the dominant two-loop correction of 0{aas) obtained 
in Refs. 



12|,[T^ and the sub-dominant corrections of 0{Gprirq 



Determining the poles of eq. (|^) corresponds to solving the equation 



Re 



^hniq'^] 



0. 



The wave function renormalization factors for the scalar Higgs bosons Hi are denoted 
as Zsi". They are the finite residues of the H and h propagators, respectively. 



^81 — ^! 



ext 
SH 



o-ext crext 

^S2 — ^Sh 



l + ReE'^^(g2) - Re 



l + ReS',,(g2)-Re 



y2 



iq' 



^q'^ -ml + T^hhiq 



^Hiq') 



-mjj + T.HH{q 




(9) 



where t'{q^) = ^t{q^). 

For a diagram with no mixing on the outgoing Higgs-boson line. Hi, the ^-matrix 
element is given by the amputated on-shell Green's function multiplied by the (^g^*)^^^. 
In the presence of mixing, i.e. h ^ H on the external line (with the scalar i being the 

^ The implementation of our calculation into HZHA is currently investigated . 
^ We do not consider here possible CT^-violating mixing between neutral Higgs bosons, which can 
occur if the MSSM Lagrangian contains complex parameters ||24|]. 
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final state particle) the respective factor reads: 

-tnniq') {Zlf)-'/' (q'-Ml 



(10) 



Therefore, in this case, the amputated Green's function is effectively multiplied by {i' 
denotes the index of the "supplementary" Higgs boson: i' = 2(1) for i = 1(2); formally 
i' = 3-i) 



Exactly analogous equations hold for the pseudoscalar constants Z^l and Z^^^. Their 
effects, however, are numerically much less important. Appropriate formulae for the 
inclusion of the Z factors into the higher-order corrected vertices can be found in the 
Appendix. 

2.3 The tteff approximation 

The inclusion of the Z factors on the external Higgs lines in the on-shell calculation 
reproduces, in the approximation of the neglected momentum dependence of the Higgs 
self-energies, the effect of using the higher-order corrected angle in an improved Born 
approximation of the cross sections (see also Ref. [^). 

The dominant contributions for the Higgs-boson self-energies (of OlGpmf) at the 
one-loop level) are obtained for = 0. Approximating the renormalized Higgs-boson 
self-energies by 

S(g2) S(0) = S (12) 
yields the Higgs-boson masses by re-diagonalizing the dressed mass matrix 

-t,H ml-t,J~^[ Ml)^ ^^^^ 

where and Mu are the corresponding higher-order-corrected Higgs-boson masses. In 
Ref. [EH] it has been shown that in the approximation with neglected external momentum 



the factors can be written as follows in terms of Aa, which is the angle required for 
the re-diagonalization in eq. (^): 

ZS'^ZrS "ff - . =+tanAa. (14) 

Z^^ = Z^l^ 5^^^^-^^ = -tanAa. (15) 

Mf^ — rrifj + L,HH 
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It is important to note that, although it is not immediately visible, both eqs. (0,^) yield 
the same angle A a ||25|| : 



Am^ + V(Am^)^ + 4S|^ 
Aa = arctan , (loj 

—2 T^hH 

where 

Am^ = (ml - ±hh) - (ml - ±hh). (17) 
The Z^^^ factor can be expressed as 

^s? ''^°2s? = cos2Aa. (18) 

l + tan^Aa ^ ^ 

The consequences for the couplings are demonstrated for Zh production in the following 
example for the ZZh vertex. The Born coupling Vzzh ~ — sin(a — /?) is changed by the 
loop corrections, in terms of the Z factors, according to 



Vzzh ~ iZ< 



1 ^ 

rcxt \ 2 



ZZh ~ {^sh) [- sin(« - /3) - ^^/l"" cos(a - (3) 

= cos Aa [— sin(Q; — /5) — (— tan Aa) cos(a — /5) ] 

= — sin(a; — /3 + Aa) 

= -sin(aefr -/5). (19) 

Analogous results hold for all Higgs vertices, including the AHi vertices. Therefore, the 
Z factors effectively shift the tree- level angle a by Aa, yielding a loop-improved angle 

aes = a + Aa (20) 

in this approximation. 

While the approximation, i.e. using an improved Born result for the cross sections 
where the tree- level angle a is replaced by cteff, incorporates the dominant one- loop and 
two-loop contributions, it is obvious from the discussion above that this approximation 
neglects many effects included in a full FD calculation. These are contributions from 
the full spectrum of the MSSM particles, the momentum dependence of the Higgs-boson 
self-energies, the gauge-boson and the fermion self-energy corrections, and in particular 
the process-specific vertex and box corrections. 



2.4 Cross sections 

In this section analytical formulae are presented for the cross sections for the on-shell 
production of the Higgs bosons e"'"e~ ZHi, e^e~ AHi including the corrections 
(i)-(iii). Box diagrams (iv) give another, more complicated, set of formfactors that make 
the expressions quite lengthy and are hence are omitted here; more details can be found 
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in Ref. [IS]. However, we include the box-diagram contributions, as described in Ref. |19 
in our numerical programs [^] and in the figures presented in this paper. 



The presented formalism for cross sections is general enough to accommodate correc- 
tions of any order to 2- and 3-point vertex functions. Beyond the one-loop level, however, 
currently only two- loop corrections to the scalar propagators have been calculated |T^,|T3|. 
Therefore, in the cross section calculations we include all possible types of one-loop cor- 
rections and the available two-loop corrections to scalar self-energies. This is well justified 
because, as discussed above, propagator corrections constitute a closed subset of the lead- 
ing 0{Gpasmf) and O^Gj^m^) contributions. Therefore, these two-loop corrections are 
of particular relevance and interest. 

The cross sections (in the CMS) for both processes and have the form: 
where A is the standard phase space factor. 



f211 



^s, m^, vr?2) = v "^^ ~'~ '^i ~'~ 2smf — Isvn^ — 2mfm2, (22) 

and ^1, A2 are defined by 

Ax + A2 cos^ ^cMs = \Y. (-^-^*) , (23) 

pol 

with M.ZS and M.ps as given below. In the following, 6'cMs denotes the scattering angle 
^CMS =^(e~,-ffj) in the CMS. The momenta of the incoming electron and positron are 
denoted as fci and ^2, respectively. The momentum of the outgoing /i, H is labeled with 
g, whereas the outgoing Z, A momentum is denoted as p, see Fig. |l]a. The matrix ele- 
ments for the Higgs-strahlung process and the associated Higgs production read (in the 
approximation of neglected box diagrams) 



u{kx)t^{p), (24) 



M 



PS 



^ ZPS \ 



Dz{s) 
'D,{s) 



v.. 



(0)Mii^7z('5) 



ZPS 



u{ki 



(25) 



For the corresponding expressions for the box contributions see Ref. . 

In the above expressions u{ki) and v{k2) are spinors of the incoming electron-positron 
pair, e^(p) is the polarization vector of the outgoing Z. cy, ca are the renormalized vector 
and axial couplings of the Z boson to an electron positron pair, at the one-loop level 

Ca = -I/Aswcw + Ca , cv = (-1 + As^^) /Aswcw + Cv\ c^y = 1 - = Mi 
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S-y(s) and I]!^^(s) denote the renormalized photon and transverse Z boson self-energies. 
Dz{s) and D^{s) are the inverse Z and photon propagators defined as 

Dz{s) = s-Ml + tl{s) , 

D^{s) = s + ±^{s) . (26) 

Finally, V denotes the eff'ective neutral Higgs-gauge-boson vertices with the one-loop 
form factors. The explicit expression for those vertices and for the matrix elements for 
Higgs-strahlung and associated Higgs production can be found in the Appendix and in 



Ref. lie 



3 Numerical results 
3.1 Parameter choice 

In the following we present numerical examples for the dependence of the neutral Higgs- 
boson couplings and cross sections on tan/5, M^, and the mixing in the scalar top sector. 
In all plots, as a typical example, the set of parameters listed in Table |^ has been used, if 
not stated differently, mj and rrih in Table |l] are the quark pole masses. M- 1 denote the 



mt 




M,- 


Ml 




Ms 


Ml 


Ma 


174.3 


4.5 


1000 


300 


200 


200 




— sLMa 



Table 1: Quark masses and SUSY parameters (in GeV) used in the numerical analysis. 



soft SUSY-breaking parameters in the scalar quark and lepton sector, respectively (in the 
following we also use the notation Msvsy = Mg) and M3 = nig denotes the gluino mass. 
The mixing in the scalar top sector, which plays a prominent role in the physics of the 
MSSM Higgs sector, is controlled by the off-diagonal term in the scalar-top mass matrix, 
mtXt = nit{At—^ cot /?), in the convention of Ref. |T^. In our analysis we have focused on 
two different values of Xt leading to two extreme values of the physical Higgs-boson mass 
Mh, as suggested by Refs. [p!3H27| -p9| . The lightest MSSM Higgs-boson mass as a function 
of Xt/MsusY has a minimum at Xf/Msusv ~ 0, denoted in the following as "no mixing" 
case. A maximum value is reached at Xf/MsusY ~ 2, denoted further as "maximal 
mixing". For the sbottom sector we assume a universal trilinear coupling. At = At. 
These values and the parameters in Table |1| are understood to be input parameters for 
the diagrammatic calculation in the on-shell renormalization scheme. 

Below we will also perform comparisons with results obtained in the framework of the 
RG improved one-loop EPA, where the input parameters are understood as MS quantities. 
To ensure consistency, in the latter case we have transformed the on-shell SUSY input 
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parameters into the corresponding MS values as discussed in Ref. pO|. The results shown 



below for the higher-order corrected Higgs-boson masses and the mixing angle within the 
RG improved one-loop EPA have been obtained with the Fortran program subhpole (based 
on Refs. MM)- 



3.2 2-loop corrections to masses and effective couplings 

The dependence of the physical neutral Higgs-boson masses on the MSSM parameters 

As an 



at the 2-loop level has been extensively discussed in the literature [l6|-|r^, |27| - 
illustration, we present in Fig. ^ the dependence of Mh on tan f] for a relatively low Ma 
value, Ma = 120 GeV. The two-loop FD result is compared with the RG improved 
one-loop EPA and also with the one-loop FD result in the no-mixing and in the maximal- 
mixing scenario in the left and in the right plot of Fig. ^, respectively. In both scenarios 
Mh shows a similar behavior: a minimum is reached around tan/? ^ 1, maximum values 
are reached for the largest tan/3 value^. In the no-mixing scenario the FD result is 
always smaller than the RG improved one- loop EPA value for Mh, with a maximum 
difference around tan /3 = 1 of up to 5 GeV. In the maximal-mixing scenario both result 
mostly agree. Note, however, that this behavior changes for larger values of Ma, where 
the maximum value of Mh obtained in the FD approach is a few GeV larger than the 
corresponding RG improved one-loop EPA value [[l^ ,|27|,^. 





taring tan/S 

Figure 2: Mh as a function of tan/3 for Ma = 120 GeV and the parameters of Table |l|. 
The no-mixing (maximal-mixing) case is shown in the left (right) plot. The solid line 
represents the two-loop FD result, the dotted line shows the RG improved one-loop EPA 
result and the dashed line shows the one-loop FD result. 



Approximate results for the cross sections have often been obtained in the literature on 



^One should keep in mind, however, that for fixed rrit — 174.3 GeV and Msusy ^ 1 TeV tan/3 
around 1 is already excluded |^,|3^ via Higgs-boson searches. 
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the basis of improved Bom results, where the effective mixing angle a^g, see eq. (]20|), and 
the higher-order corrected Higgs-boson masses and have been evaluated within 
the RG improved one-loop EPA. We will in the following compare our FD result for the 
cross sections with this approximation, to which we will refer as "RG aes approximation" . 
As mentioned above, the results in the RG aes approximation have been obtained using 
the program subhpole (based on Refs. [0,0,0). 

In order to disentangle the effect of different contributions in this comparison, we 
first discuss the results for the effective Z-Higgs-boson couplings in the two approaches. 
In Sect. p^.3| we have shown that the contribution of the wave function renormalization 
factors of the Higgs bosons is given by the effective mixing angle ctefr evaluated in the 
FD approach, if the momentum dependence of the Higgs-boson self-energies is neglected. 
However, in the actual cross-section calculation in the FD approach the momentum de- 
pendence is included at the available (currently one-loop) level. Therefore, for a better 
comparison of the quantities actually entering the cross section calculation in the two 
approaches, we formally define aes in the FD approach as (in analogy to eq. (p!5|)) 

a,\ = -arctanZsT, (27) 

where the Z^f^ is given by the exact expression (p!TD (with Higgs self-energies calculated 
at = Ml), not by the = approximation of eq. (|l5l)p|- 

Using this definition, in Fig. |^ the dependence of sin(aeff — /3) on tan /? and is shown 
in both approaches, for M4 = 100, 120, 150 GeV and the no-mixing scenario, Xf = 0. 
Since Mh is a derived quantity and not an input parameter in our approach, the parameter 
that is actually varied in the plots shown as function of Mh is tan /3. For simplicity, i.e. 
in order to avoid a non-functional behavior, in all plots shown as function of Mh in this 
paper we restrict the tan j3 region to tan (3 > 1.6 (as mentioned above, for rrit = 174.3 GeV 
and MsusY ^ 1 TeV tan P values around 1 are already excluded [^|31[] via Higgs-boson 
searches.) 

As can be seen in Fig. ^ the agreement between the FD two-loop result for a^g and 
the result within the RG improved one-loop EPA is in general remarkably good (while 
large deviations can appear compared to the FD one-loop result, see e.g. the middle plots 
in Fig. For most of the tan/3 range the FD two- loop result and the result within the 
RG improved one- loop EPA differ by not more than 5%, larger deviations can be observed 
only for Ma = 100 GeV and tan/3 > 10 (left upper plot in Fig. H), where sin(aeff — P) 
itself is small. Even in this case the sin(aeff — P) values in both approaches agree very well 
with each other when expressed in terms of the physical Higgs mass Mh- For the maximal- 
mixing case, Xt/MsusY = 2, the differences between the effective couplings obtained in 
both methods are even smaller. 

It should be noted that the behavior of sin(aeff — /3) in the limit of large tan P is quite 
different for small and large pseudoscalar masses. This behavior changes for M4 between 
100 and 150 GeV; the actual value depends on the stop mixing parameter Xt (see also 
Fig. 9 in Ref. pSl.) 



^ It should be noted that defined as = arctan2^™|f shghtly differs from a^g. 
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Figure 3: sin(Q;efr — /3) as a function of tan/3 (left plots) and Mh (right plots) for Xt = 
and the parameters of Table |^. The solid line represents the two-loop FD result (see 
eq. (p7D ), the dotted line shows the RG improved one-loop EPA result, and the dashed 
line shows the one-loop FD result. 
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3.3 Results for the cross sections 



Differences in the Higgs-production cross sections between our FD result (containing the 
complete one-loop result and the dominant two-loop corrections) and the RG acs ap- 
proximation have a two-fold origin: the different predictions for the values of Mh and 
Oes, which we compared in the previous section, and the additional contributions con- 
tained in the FD result (i.e. the one-loop 3- and 4-point vertex functions and the 2-point 
contributions that are not contained in Oes)- 

In Fig. ^ we present the cross sections for the two production channels for a LEP2 
energy of ^/s = 206 GeV, in the no-mixing and the maximal-mixing scenario, as a function 
of tan/3. Fig. ^ shows the same results as a function of M^. At LEP2 energies, the box 
diagram contributions are small, of the order of 2-3% and do not modify the cross 

section behavior in a significant way. 

As can be seen from both figures, the cross sections for the Higgs-strahlung process 
corresponding to the FD result and the RG approximation are close to each other, 
with differences of the order of a few per cent. The only exception occurs at f» and 
large tan (3. In this case the difference can amount up to 50% when the cross sections are 
expressed as a function of tan/3, but is considerably reduced when they are calculated 
in terms of Mh- The two-loop result in the FD calculation is always above the one- loop 
result; the difference can be sizable, but is then mostly due to the kinematical effect that 
the values for at one-loop are much larger than at the two-loop level. This effect is 
especially pronounced at the kinematical endpoints, i.e. for large tan p. 

Associated Ah production is of interest at LEP2 energies only for sufficiently low 
Ma_ < 120 GeV, otherwise it is kinematically forbidden. Therefore we restricted our plots 
for a^h in Figs. ^ and |^ to Ma = 100 GeV. The difference between the FD result and 
the RG tteff approximation is larger than for azh and may reach about 20% for not too 
large tan f3 values. The result of the two-loop FD calculation is mostly above the one-loop 
result, but the differences are much smaller than for the Higgs-strahlung process. 

For both processes, the differences between the FD result and the RG acs approxi- 
mation are not particularly pronounced, however visibly larger than those induced by the 
modifications of ttefr only. For LEP2 energies they can be attributed mostly to the vertex 
corrections to the ZZh and ZhA couplings . 



The dependence on the t mixing is depicted in Fig. |^. The main effect on azh, (^Ah is 
kinematical and arises from the change of with Xt/ Msusy ■ This effect leads to a visible 
(additional) asymmetry in the two-loop FD results, whereas the RG Oefr approximation 
and the one- loop FD result show a weaker asymmetry in Xt/ Msusy- As can be seen 
from the figure, the differences between the methods are much larger (typically more then 
10%) when associated production is considered, with the above-mentioned "kinematical" 
asymmetry further increased by the inclusion of the vertex corrections. 

In Fig. 1^ the cross sections for the Higgs-strahlung process and the associated produc- 
tion are shown for a typical Linear Collider energy, ^/s = 500 GeV 0, in the maximal- 
mixing scenario (in the no-mixing scenario similar results have been obtained). In addition 
to the previous plots we also show the result for the two-loop FD calculation where the 
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Figure 4: cr^/i and a^/i at ^/s = 206 GeV as a function of tan/3 for two values of 
and the parameters of Table |l]. The upper (lower) row contains the result for the no- 
(maximal)-mixing scenario. The solid line represents the two- loop FD result, the dotted 
line shows the result for the RG aes approximation and the dashed line shows the one-loop 
FD result. 



box contributions have not been included, in order to point out their relative importance 
for high-energy collisions. For ^/s = 500 GeV the differences between the FD result and 
the RG OefT approximation are larger than in the low-energy scenario, an effect that is 
more pronounced for the higher value of M4. For Ma = 200 GeV, typically they are 
of the order of 10-15% for azh and even up to 25% for a ah (for ^ 90 GeV). The 
difference between the two-loop and one-loop FD result can be sizable. The two-loop 
result for azh is in general larger than the one- loop value, again increasing with Ma, 
where for Ma = 200 GeV the difference can amount up to 15%. aAh, on the other hand, 
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Figure 5: azh and a Ah at \/s = 206 GeV as a function of for two values of Ma and the 
parameters of Table |I]. The upper (lower) row contains the result for the no- (maximal)- 
mixing scenario. The solid line represents the two-loop FD result, the dotted line shows 
the result for the RG aefr approximation and the dashed line shows the one-loop FD 
result. 



is decreased at the two-loop level for Ma = 200 GeV and the difference may be sizable. 

The box contributions become more important for higher ^/s and change the total 
cross section by 5-10%. This result remains unchanged even if sleptons are significantly 
heavier than Mj = 300 GeV used in our numerical analysis, as the dominant contributions 
to box diagrams are given by W and Higgs boson exchanges |jl8|,|l9|, which do not depend 
on Ml. Also, one should recall that box contributions lead to an angular distribution 
of the final-state particles different from the effective Born approximation and thus give 
much larger corrections to the differential rather than to the total cross section, at least 
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Figure 6: azh and a Ah as a function of Xt/MsusY for = 206 GeV, tan/? = 5 and the 
parameters of Table ^ The sohd hne represents the two- loop FD result, the dotted line 
shows the RG Oefi approximation and the dashed line shows the one- loop FD result. 



for some range of the scattering angle. Therefore, the box diagrams have a significant 
effect at Linear Collider energies and thus have to be included. The same conclusions can 
be drawn for ^/s = 500 GeV in the no-mixing scenario, which we do not show here. The 
differences between the FD result and the RG aes approximation are only slightly smaller 
than in the maximal-mixing case. 

In Fig. p the results for azh and a Ah are shown as a function of ^/s in the no-mixing 
scenario for Ma = 100, 200 GeV. Besides the obvious kinematical drop-off of the cross 
sections, one can observe that the relative differences between the FD two-loop result and 
the RG Oeff approximation grow with ^/s. The differences remain almost constant or even 
increase slowly in absolute terms, whereas the full cross sections decrease. aAt becomes 
very small for large Ma, as can be seen in more detail in Fig. ^. There we show the 
dependence of azt and (Tah on Ma in the no-mixing scenario for tan (3 = 5 and tan (3 = 50. 
For azh, the A boson decouples quickly; the dependence on Ma becomes very weak for 
Ma ^ 250 GeV, when azh is already practically constant (compare e.g. Ref. [Il3[). In the 
same limit, a Ah goes quickly to zero due to suppression of the effective ZhA coupling, 
which is ~ cos(aefr — (3); also the kinematical suppression plays a role, but this becomes 
significant only for sufficiently large Ma, Ma > 350 GeV. For large tan/? the decoupling 
of Ma is even more rapid. The differences between the FD two-loop result and the RG 
tteff approximation for the Higgs-strahlung cross section tend also to a constant, but they 
increase with Ma for the associated production. The latter can be explained by the 
growing relative importance of 3- and 4-point vertex function contributions compared to 
the strongly suppressed Born-like diagrams. As can be seen from Fig. |^, for tan (3 = 50 
and Ma > 300 GeV the FD two-loop result is almost an order of magnitude larger than 
the result of the RG a^s approximation, and starts to saturate. This can be attributed 
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Figure 7: azh and aAh as a function of at ^/s = 500 GeV, for two values of Ma, 
Xt/MsusY = 2, and the parameters of Table |l]. The solid (dot-dashed) line represents 
the two-loop FD result including (excluding) box contributions, the dotted line shows the 
RG Q!efr approximation and the dashed line shows the one- loop FD result. 



to the fact that the (non-decoupling) vertex and box contributions begin to dominate 
the cross section value. However, such a situation occurs only for very small a Ah values, 
aAh ~ 10~^ fb, below the expected experimental LC sensitivities. 

4 Conclusions 

Using the Feynman-diagrammatic approach we have calculated the production cross sec- 
tions for e^e~ — hZ, hA. The Higgs-boson propagator corrections have been evaluated, 
including besides the full one-loop result also the dominant and subdominant two-loop 
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Figure 8: azh and a Ah as a function of i/i for tan/3 = 5, Xt = 0, and the parameters of 
Table ^ The sohd (dot-dashed) hne represents the two-loop FD result including (exclud- 
ing) box contributions, the dotted line shows the RG acfr approximation and the dashed 
line shows the one- loop FD result. 



corrections. In addition we have included also the full set of one-loop vertex and box 
corrections. 

We have also investigated an improved Born approximation based on the effective 
mixing angle, Oefr, in the neutral CP-even Higgs sector. We have shown analytically 
that this approximation corresponds to taking into account the Higgs-boson propagator 
corrections with neglected external momentum at one-loop and two-loop order. 

While the previous results available in the literature neglected the large two-loop cor- 
rections in the case of the Feynman-diagrammatic approach and were restricted to an 
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Figure 9: azh and a Ah as a function of at ^/s = 500 GeV, shown for tan /5 = 5 and 
tan j3 = 50, Xf = 0, and the parameters of Table |I]. The sohd (dot-dashed) hne represents 
the two-loop FD result including (excluding) box contributions, the dotted line shows the 
RG tteff approximation and the dashed line shows the one- loop FD result. 



improved Born approximation in the case of the renormalization-group improved one- 
loop Effective Potential approach, our new result combines the dominant two-loop correc- 
tions with the complete Feynman-diagrammatic one-loop result and thus represents the 
presently most precise prediction for the production cross sections for e^e~ — >■ hZ, hA. 

Specifically we have numerically analyzed the effect of the two-loop contributions 
incorporated in our result, which turned out to be sizable. We have furthermore compared 
our full Feynman-diagrammatic results for azh and aAh with the approximation where 
tteff and Mh are evaluated within the renormalization-group improved one-loop Effective 
Potential approach (RG Ocs approximation). For LEP2 energies, ^/s = 206 GeV, the 
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difference between the Feynman-diagrammatic result and the RG acfr approximation stays 
mostly at the per cent level in the parameter space allowed by the LEP2 exclusion limits 
for Mh. At energies reachable at an e+e" linear collider, e.g. ^/i = 500 GeV, the difference 
between the Feynman-diagrammatic result and the RG a^s approximation can reach 10- 
15% for (Tzh or even 25% for aAh- The box contributions play an important role for 
high ^/s and can amount up to 10%. Therefore in a precision analysis for a high energy 
e+e" collider the two-loop propagator corrections as well as the complete diagrammatic 
one-loop contributions should be included]^. 
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Appendix Explicit expressions for the cross sections 

Performing the sum over polarization states in eq. ( ^3]) one gets for the associated scalar 
and pseudoscalar production: 



IPS 



- A: 



2PS 



''PS 



+ 



'PS 



(A.l) 



where 
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2V: 



ZPS 



(A.2) 



and A as given in eq. 
are more complicated: 



A 



2ZS 



The corresponding expressions for the Higgs-strahlung process 
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The quantities F and G denote effective Higgs-gauge-boson vertex form factors. They 
can be defined as follows. At the tree level the relevant Higgs-boson vertices read (the 
assignment of momenta is given in Sect. 
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where V^zs ^^^1 V^p^ can be written in matrix form as (index j numerates CP-odd Higgs 
bosons: A = Pi, G = P2): 



v.. 



{0)i 



zzs 



eMz ( cos(a - (3) 



(A. 



v.. 



(0)ii 



sin(a — [3) cos(a; — [3) 



ZPS 



(A.9) 



2swCw \ cos(a -13) - sin(a; - [3) j ' 
After including one-loop vertex corrections the renormalized vertices have the form: 

(A. 10) 



^ZPS - P [ ^ZPS + 



^7'^^ {V^% + Fl) + p^p'F^ + q^^q^Fl + p^g^F] + q^p^F^^ , (A.ll 



where F^'s are the renormalized vertex form factors. 

Next, we define the respective quantities with tilde which are obtained by inclusion 
of the one-loop corrections on the external lines. For instance, the complete vertices read 
(in the following: i' = 3 — i, j' = 3 — j.): 
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(A.14) 



The vertices Vzps, ^zzs can be decomposed into V^p^, Pp, -F5 etc. (as in eq. ( |A.l(j| , |701] )) 
where e.g. F{ = (Zf'^ZIf )^/^ (F{ + Z^I^F() etc. 

To include photon exchange in the s-channel one needs to know also the renormalized 
vertices V^ps and V^zs (which vanish at the tree level): 



V^pg — P ^ p q (J 



5' 



(A.15) 



gf^'Gl + p^p^Gl + q^q^Gl + p^q^Gl + q^p'G' 



(A.16) 



The vertices V^ps, V^zs are defined similarly as shown in eqs. ( |A.12| , |A.13| ). 

The explicit expressions for the one- loop corrections to the vertices VzpS) Vzzs, V'yPS 
and Vyzs can be found in Ref. |TB . 
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